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Electromagnetic Hadronic Form-Factors 

Robert G. Edwards^, for the Lattice Hadron Physics Collaboration 

'^Thomas Jefferson National Accelerator Facility, Newport News, VA 23606, USA 

We present a calculation of the nucleon electromagnetic form-factors as well as the pion and rho to pion 
transition form-factors in a hybrid calculation with domain wall valence quarks and improved staggered (Asqtad) 
sea quarks. 
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The determination of light hadron physics 
properties is of considerable interest to experi- 
mental labs such as Jefferson Lab. Many quanti- 
ties, such as hadron form-factors and generalized 
structure functions, are not well understood the- 
oretically. The available and forthcoming high 
precision experimental data presents consider- 
able challenges and many opportunities for corre- 
sponding lattice calculations. Including the con- 
tribution of light dynamical quarks is an impor- 
tant goal of this work, but the cost of generating 
these gauge ensembles is large. We have adopted 
a so-called "hybrid" scheme where staggered sea 
quarks are used (Asqtad action ^) and domain 
wall valence quarks P|. While unitarity is bro- 
ken at finite lattice spacing, it is recovered in the 
continuum limit when the sea and valence quark 
masses are properly tuned. 

In this contribution, we study the efficacy 
of an uncommon method of calculating three- 
point functions that avoids sequential source tech- 
niques. We calculate various electro-magnetic 
form-factors without the cost of a new sequential 
inversion for each new set of observables. The 
method is particularly suitable for valence quark 
actions amenable to a multi-mass solver, such as 
the Overlap quark action. The conclusion is that 
the method appears viable, but very light quark 
mass tests are needed (and on-going). 

A typical sequential source method for comput- 
ing a matrix element such as {B{pf )\V{q)\A{pi)) 
for an initial state A with three-momenta pi and 
final state B at momenta pj with a two-quark 
insertion V involves either a sequential inversion 
through the insertion or the sink. To properly 
extract the matrix element for a set of momenta 



q = Pf — Pi, usually the sequential source for the 
sequential insertion is held at a fixed momenta. 
There are two common techniques. Sequential 
inversion through insertion: benefits are that one 
can vary the source and sink fields but the inser- 
tion momenta and operator are fixed. Sequential 
inversion through sink: benefits are one can vary 
the insertion operator and momenta, but the sink 
operator and momenta are fixed. In addition, for 
baryons the spin projection matrix between the 
source and sink must be fixed necessitating dif- 
ferent sequential inversions to extract electric and 
magnetic quantities. The common problem is one 
vertex must have a definite momentum. 

As an alternative, one can instead make a sink 
(or source) quark propagator (as opposed to the 
state) have a definite momentum. Putting the 
sink (or source) quarks at definite momentum im- 
plies the hadron state has a fixed momenta via the 
translation/momentum operator on the lattice. 
The scheme then is to build the desired hadron 
state propagating from the source and sink where 
the quark propagators used for one of these states 
has a fixed momenta. Thus, one avoids sequential 
inversions computing {B{pf)\V{q)\A{pi)). For 
example, fixing the sink momenta to a fixed pf, 
one can extract the matrix element at all q by 
vary pi according to momentum conservation. 

To fix a quark's momenta necessarily implies 
a wall source thus requiring the need for gauge 
fixing. Additional tricks can be used to improve 
statistics like using charge conjugation and time 
reversal (CT) in (anti-)periodic boundary con- 
ditions. This method also works for Dirichlet 
boundary conditions where one maintains equal 
source and sink separation from the Dirichlet 
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Figure 1. Pion electromagnetic form- factor com- 
puted by Eq.nand imposing the lattice dispersion 
relation. 



Figure 2. p n transition form-factor. The 
VMD line has a 0.4 normalization. 



wall. E.g., a wall sink hadron because a wall 
source after CT. If applicable, multi-mass inver- 
sion techniques can be applied both to the source 
and the sink propagator calculations. 

For our calculation of hadronic form factors 
with unqucnched gauge configurations, we per- 
formed a hybrid calculation using MILC Nf — 
2 + 1 configurations in 20'^x64 volumes, gener- 
ated with staggered asqtad sea ciuarksT, and 
domain wall valence quarks with domain wall 
height mo = 1.7 and extent Ls = 16 of the 
extra dimension The MILC configurations 
were HYP blocked |Sj before valence propaga- 
tors were computed, otherwise the residual chiral 
symmetry breaking would have been unaccept- 
ably large. Dirichlct boundary conditions were 
imposed 32 time slices apart. The lattice spac- 
ing scale is set at a = 0.125fm. The Asqtad set 
has m„ = 0.01 and uis = 0.05, P = 6.76. The 
valence pion and rho mass are tUt^ = 758(5)MeV 
and TOp = 1060(14)MeV, resp. 

The technique for extracting form factors uses 
the ratio method of correlation functions: 

(1) 



where A, B, C are generic smearing labels, L is lo- 
cal and fi is the direction of the local current, the 
initial and final states are Xa and Yp. Similarly 
for R^^ = R^^ where Y ^ X . Note, momenta 
and smearing labels are not interchanged. In the 
case oiY ^ X, this ratio cancels all wave- function 
factors and exponentials. For the transition case, 
the combination {R^-^ R^^y^^ cancels all expo- 
nentials and wave- function factors taking into the 
effects of smearing on the wave-functions. One 
extracts the form-factors up to some kinematic 
factors using properties of the spin projection ma- 
trix between Xa and Yg. In all the subsequent 
work, we fixed the wall sink to have pf = 0. The 
source quarks are APE smeared. 

The first test of the wall-sink method involves 
the pion form factor. The result shown in Fig. 
compares favorably with a sequential-sink tech- 
nique 0] at the one mass in common correspond- 
ing to arrival = 0.081. The obtained value of Zy 
for our local current also agrees well. 

Electro-disintegration of the deuteron has been 
intensively studied experimentally. While isovec- 
tor meson exchange currents have been identified 
in such systems, the role of isoscalar meson ex- 
change is not clear. Here, we provide a lattice 
measurement of the p — > tt transition form fac- 
tor which is an important step in an analysis of 
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DWF Proton Magnetic Form-Factor - Wall Sink 
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Figure 3. Proton electric and magnetic form fac- 
tors. Also shown are results from QCDSF [Jj. 



exchange mechanisms. The matrix element com- 
puted is 

iriTj + mp 

The form factor V{Q'^ ) shown in Fig.[21is in rough 
agreement with a phcnomenological calculation 
in Ref. j^j- The coupling constant normalization 
taken from 1^(0) is in rough agreement with ex- 
periment Ref. |6I5| . Higher are needed. 

The extraction of the proton and neutron form 
factors uses Eq. ^ together with a spin projec- 
tion matrix Tap between source and sink nucleon 
states Na and Np. We can increase statistics on 
the electric and magnetic form-factors by averag- 
ing over all the nucleon spin polarizations. The 
results for the proton and neutron are shown in 
Figs. 13 and ^ Compared with quenched clover 
results from QCDSF '7 at roughly the same pion 
mass (/? = 6.0, n = 0.1338), the error bars stay 
somewhat constant at higher . Reasonable 
agreement is seen with experiment and Ref. 
In addition, a signal is seen for the computation- 
ally difficult neutron electric form-factor. 

This work was supported by the U.S. De- 
partment of Energy under contract DE-AC05- 
84ER40150. Computations were performed on 
the Pentium IV clusters at JLab. 
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Figure 4. Neutron magnetic form-factor. 



